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Contribution of order a 2 {Zct) & \v?{Za)~ 2 to the ground-state Lamb shift 
in hydrogen induced by the loop-after-loop diagram is evaluated analytically. 
An additional contribution of this order is found compared to the previous cal- 
culation by Karshenboim [JETP 76, 541 (1993)]. As a result, an agreement is 
achieved for this correction between different numerical and analytical meth- 
ods. 

In the present work a part of the two-loop self-energy correction to the Lamb shift in 
hydrogen is investigated, which is induced by the irreducible part of the diagram in Fig. 1(a). 
It is referred to as the loop- after-loop contribution. The second-order self-energy correction 
is important for the comparison of the theoretical prediction with the experimental results 
for the Lamb shift in hydrogen Jl| and He + @ and, therefore, it influences the charge radius 
of the proton which can be extracted from the Lamb shift in hydrogen [[| . 

The loop-after-loop contribution has been the subject of a recent debate in the litera- 
ture. Analytic calculations of its Za-expansion coefficients have been carried out by Eides 
and co-workers |4[] and Pachucki [|[ in order a 2 (Za) 5 and by Karshenboim ||[] in order 
a 2 (Za) 6 ln 3 (Za)~ 2 . A direct numerical calculation of this correction to all orders in Za 
in the low-Z region was reported by Mallampalli and Sapirstein ||. A fit to the data in 
Ref. H is consistent with the analytical result in order a 2 (Za) 5 but it is in significant dis- 

1 The same result was re-derived recently in Ref. |?J 
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agreement with Karshenboim's calculation of the leading logarithmic contribution. As a 
consequence, a question was raised in Ref. about the breakdown of the Za expansion for 
the two-loop self energy even for hydrogen. The subsequent calculation by Goidenko et al. 
H, also nonperturbative in Za, shows to be compatible with the analytical results. How- 
ever, my recent nonperturbative (in Za) calculation [|TtJ exhibits a good agreement with the 
results of Ref. f| and yields a corresponding logarithmic contribution which is roughly three 
times larger than Karshenboim's result. On the contrary, a recent evaluation by Manohar 
and Stewart [|ll]] supports Karshenboim's value of the leading logarithmic correction for the 



total two-loop self energy. 

In Refs. ft is argued that the leading logarithmic correction for the two-loop self en- 
ergy is induced only by the diagram shown in Fig. 2(b), if working in the Fried- Yennie gauge. 



The numerical calculation of Ref. |T0[ confirms their value of the corresponding contribution 
from the diagram in Fig. 2(b). However, it shows also an additional logarithm cubed with 
a coefficient = —0.652(30) originating from the diagram presented in Fig. 2(a). In the 
present investigation, I derive the leading logarithmic contribution induced by this diagram 
in the case of the ground state of hydrogenlike ions. The resulting coefficient is found to be 

= —2/3, in agreement with my previous numerical result. 

The contribution of the loop-after-loop diagram [Fig. 1(a)] can be written as 

A£l a , = I dpdkdk!dp'$ a (p)V R (j>, k)G r a cd (k, k')X«(k', p / )Va(p') , (1) 

where ip a is the ground-state wave function, E# is the renormalized self-energy operator, 
G„ ed is the reduced Coulomb Green function, and the time component of all 4-vectors equals 
to the ground-state energy e a . This expression is not completely gauge invariant, but it 
can be shown to be gauge independent over covariant gauges. In the present investigation 



2 Relativistic units are used in this Letter (h = c = m = 1). I use roman style (p) for 4-vectors, 
bold face (p) for 3-vectors and italic style (p) for scalars. 4-vectors have the form p = (po, p). The 
length of the 3- vector p is denoted by p. I also use the notations ]f> = p M 7 At , p = p/|p| . 
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the consideration will be given in the Fried- Yennie gauge. Because of lack of a convenient 
representation of the Coulomb Green function, the standard way to evaluate Eq. ([[J) is to 
expand the inner electron propagators in terms of an interaction with the external nuclear 
field. Some of the terms of this expansion are shown in Figs. 2[(a),(b)]. 

The leading contribution to AE\ a \ is of order a 2 (Za) 5 and originates from the diagram 
in Fig. 2(a). It was evaluated in Refs. after proving that only relativistic values of 
the intermediate momentum are responsible for this correction and that both external mo- 
menta can be put on the mass shell with vanishing space components. The next-to-leading 
contribution contains a logarithm cubed. It was derived in Refs. ||[7|] by using an effective- 
potential approach, where the self-energy operator is replaced by a potential on the mass 
shell, 

/ , s t t 8a(Za) , , 
S R (p, k) -> Vse = In k . (2) 

The effective quasilocal potential Vse is defined in a way to generate the leading contribu- 
tion to the one-loop Lamb shift. Obviously, this potential leads to a squared logarithmic 
contribution in Eq. ([!]), after the integration over the intermediate momentum is carried 
out within the logarithmic region Za < < 1. One can gain one more logarithm if In 2 k is 
integrated with a factor 1/k, 

- 1 \n 2 k 



dk^r 11 ~ In 3 Za . (3) 

Za k 

In Refs. the factor 1/k arises from the second term of the potential expansion of G™ d 
in Eq. ([I]). This is the reason why in these considerations Fig. 2(b) yields a logarithm cubed 
and Fig. 2(a) does not. As we will see below, the treatment of Refs. ||[7| is incomplete. A 
more accurate investigation shows that an analogous factor 1/k can be obtained by taking 
into account the momentum distribution of the external wave functions in the diagram in 
Fig. 2(a). 

Let us now derive the next-to-leading contribution to the diagram in Fig. 2(a). While its 
leading term can be evaluated assuming the external momenta being put on the mass shell 



with vanishing space components, we intend to keep the first-order terms in the expansion 
over the space components. The binding energy is still neglected; one can show that it does 
not contribute to the order of interest. We write the contribution of the diagram in Fig. 2(a) 
as 

AE2a = (2^ / ^W^I(P)^(P - k)/C(p, k, p')V c (k - p')^(p') , (4) 

where Vc is the Coulomb potential Vc(q) = —AnZa/q 2 , and K, is the kernel 

/C(p,k,p') = 7 oA^(p,k)i^-A^(k,p') , (5) 

with being the renormalized vertex operator in the Fried- Yennie gauge. For our purposes 
it is sufficient to keep only the lowest-order term of the Za expansion of the Dirac wave 
function 

I ( ) = ( 9aijp)Xnam a {v) \ _ N ( X^nia (p) \ 

Wa{P ' \fa(p)X- Ka mM ) ~ (P 2 + (ZaW \(-p/2) X - Ka mM J ' 1 ' 

where N = (2vr) 3 / 2 (32(Z«) 5 /vr) 1/2 . In Eq. (|) and in what follows, the sign of the approxi- 
mate equality (~) is used when terms irrelevant for the contribution of interest are dropped. 
Assuming that the external momentum p is close to the mass shell, we can expand the vertex 
operator over the space components of p 

Agf ( P) k ) « A (°) ( S , k ) + (7 . p ) A « (s, k) + (p • k) A< 2) (s, k) + p 2 A< 3) (s, k) , (7) 

where s is a 4- vector with zero space components, s = (1,0). The parts of AE2 a induced by 
the four terms in the righthand side of Eq. (^) are denoted by AEi, i = . . . 3. 

The correction AE can be written by introducing an effective wave function ip as 

AE = J ^-^t (k)/c(g) k; s)y , o(k) j (8) 

where 

dp N ( X Ka m a (k) ^ 



Vo(k)=/^^(k- P )^(p) - 2k2 



Za k . - . 



(9) 



Eq. (U) contains a contribution of the previous order. To eliminate it, one should subtract 
the same expression with the Dirac bispinors replaced by the Schrodinger wave functions. 
This means that we are interested only in the part of Eq. @ which contains a product of 
the upper and lower components of <£>o- Considering the matrix structure of the kernel /C, 
one can see that only the part of K, proportional to 70(7 • k) provides a nonzero contribution 
to the order of interest. Therefore, we have 

^p^)^- ( rj )2 7o( p k) ^w- (io) 

One can show that J-'o(k) « — (112/3) A: 2 In 2 k 2 for small k. Keeping in mind that logarithms 
normally arise from integrals like Eq. (^), we can cut off the integration over k in Eq. (^) 
by an arbitrary constant A ~ 1 and use the low-energy asymptotic form for ^(k), 

a 2 {Zaf\Yl M „ In 2 k 2 k 

AE « hr^^T / dk——SiTCtg— . 11 

2n 6 3 Jo k Za 

If Za is set to zero under the integral, it would logarithmically diverge for small values of 
k. Therefore, Za can be regarded as an effective cut-off parameter in this case. We can cut 
off the fc-integration by Za and then set Za to zero under the integral 

f A „ln 2 fc 2 k 71 f A , In 2 A; 2 ir , , „ 

J = / dk^— arctg— w - / dk^— w — In 3 (Za " 2 . 12 
io k & Za 2 Jza k 12 v 7 

A rigorous evaluation of the integral J yields the same value of the cubed logarithmic 
contribution. So, we have for the correction AEq 

AE « (-^) ^(Za) 6 ln 3 (Z«)^ 2 . (13) 

The evaluation of the remaining corrections is carried out in a similar way. By a direct 
calculation one can show that AEi w 0. The correction AE 2 is written as 

rfk 



AE- 
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J^M^ 2 (k)7oA ° 2) (S ' k) (K + 1 } A ° 0) (k ' S) Vo (k) 



+^ (k) 7 oAi 0) (s,k)(K+ l)A( 2) (k,s)^ 2 (k)} , (14) 
with a wave function (p 2 



w(k) = / ^5(p-k)li(k-pW.(p) 



V a )a r ctgA(w».(k)\ (15) 



2A; 2V y & Za V 
Keeping only terms proportional to the product of the upper components, we have 

AE 2 » 4^- / rfA;^^ arctg— , 16 

where the function T 2 {k) can be found to have an asymptotic form T 2 {k) ks — (64/9)A; 3 In 2 k 2 
for small fc. Immediately, we have 



AE 2 « | ^(Zc*) 6 J « g ^(Zc*) 6 ln 3 (zTa)- 2 . (17) 



The correction Ai? 3 is evaluated in the same way as AE 2 by introducing a new wave function 

dp 



(27T) 

- --(Za) arctg— ^ Q j. (18) 

The derivation yields 

AE 3 ^^-^(Zaf\n 3 (Za)- 2 . (19) 

Adding all evaluated corrections, we have for the total cubed logarithmic contribution from 
the diagram in Fig. 2(a) 

3 / 2\ n 2 

AE 2a ^J2 AE * « ho ^(^«) 6 ln 3 (zTa)- 2 . (20) 

i=0 \ O/ 7T 

The results are commonly represented in the form of a double expansion over (Za) and 
ln(zTa)- 2 

AE = (f) 2 XX (^)'w(^)" 2 • (21) 

In this notation the leading logarithmic contribution induced by the diagram in Fig. 2(a) is 
found to be = —2/3, in a good agreement with a 63 = —0.652(30) obtained in my previous 
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calculation fllQ| . This result, together with the corresponding contribution from the diagram 
in Fig. 2(b) evaluated in Ref. ||, yields the total correction of order a 2 (Za) 6 \n s (Za)~ 2 for 
the loop-after-loop diagram [Fig. 1(a)], a 63 = -2/3 - 8/27 = -26/27 w -0.963. This value 
is in good agreement with the numerical result of Mallampalli and Sapirstein = —0.9 || 
and with my numerical calculation = —1.01(8) fllPf . However, it strongly disagrees with 



the evaluation by Goidenko and co-workers |§. A possible reason for this discrepancy may 
be a non-covariant renormalization procedure used in that work, which is known to provide 
spurious contributions in some cases. This topic is discussed in more detail in Ref. |0[. 

In the present work, only the irreducible part of the diagram in Fig. 1(a) is investigated. 
The question if the remaining two diagrams of the two-loop self energy provide a cubed 
logarithmic contribution, remains still unclear. An indication that it could be the case are 



the results of Pachucki |12] and Manohar and Stewart [[Tl|]. The authors calculated the 
leading logarithmic contribution to the total two-loop self energy by working in Coulomb 
gauge [|X2j] and within the renormalization-group approach ||11|| . Both evaluations agree with 



Karshenboim's result. Keeping in mind that we derived the additional contribution from the 
diagram Fig. 1(a), one should deduce that the remaining diagrams Figs. l[(b),(c)] provide 
the same contribution with the opposite sign, if working in a covariant gauge. 

In summary, I derived an additional cubed logarithmic contribution to the loop-after-loop 
correction for the ground state of hydrogenlike ions. As a result, an agreement between the 
results of numerical and analytical methods is achieved. Therefore, there are no grounds at 
present to question the validity of the Za expansion for hydrogen. However, the convergence 



of this expansion for the two-loop self energy is remarkably slow. As was shown in Ref. |T0 
the first three expansion terms for the loop-after-loop correction cover only about 50% of 
the total result for hydrogen. Thus, it is highly desirable to carry out exact calculations of 
the remaining part of the two-loop self energy in the low-Z region. 
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FIG. 1. One-electron self-energy Feynman diagrams of second order in a. 
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FIG. 2. Two terms of the potential expansion of the diagram in Fig. 1(a). A double line denotes 
an electron in the field of the nucleus. A single line indicates a free electron. A dashed line denotes 
a Coulomb interaction with the nucleus. 



9 



